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Abstract. We investigate the possibility of measuring the Hubble constant, the fractional energy density com- 
ponents and the equation of state parameter of the "dark energy" using lensed multiple images of high-redshift 
supernovae. With future instruments, such as the SNAP and NGST satellites, it will become possible to ob- 
serve several hundred lensed core-collapse supernovae with multiple images. Accurate measurements of the image 
separation, flux-ratio, time-delay and lensing foreground galaxy will provide complementary information to the 
cosmological tests based on, e.g., the magnitude-redshift relation of Type la supernovae, especially with regards 
to the Hubble parameter that could be measured with a statistical uncertainty at the one percent level. Assuming 
a fiat universe, the statistical uncertainty on the mass density is found to be of^T JJ0.05. However, systematic 
effects from the uncertainty of the lens modeling are likely to dominate. E.g., if the lensing galaxies are extremely 
compact but are (erroneously) modeled as singular isothermal spheres, the mass density is biased by ct'q 3 ~ 0.1. 
We argue that wide-field near-IR instruments such as the one proposed for the SNAP mission are critical for 
collecting large statistics of lensed supernovae. 
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1. Introduction 

Gravitational lensing of high-z objects has been used in 
the past as a tool for deriving cosmological parameters 
with mixed success. While the fraction of quasars with 
multiple images and the distribution of image separa- 
tions may be used to probe the vacuum energy density, 
Q A ([Turner 1990|; [Schneider et al. 1992|), the method suf- 



fers from severe systematic uncertainties related to the 
lens modeling and the results are as of yet inconclusive 



(Fukugita et al. 1992 


Maoz & Rix 1993|; Kochancck 1996: 


Park & Gott 1997|; ( 


]hiba & Yoshii 1999|). On the other 


hand, the method proposed by 


Rcfsdal (1964 1966|) using 



time delay measurements of multiple imaged quasars to 
constrain H has provided results that are in good agree- 



ment with independent techniques, (see e.g. Koopmans 



& Fassnacht 1999; Brown 200C). Thus, distances derived 



from geometrical measurements of lensed high-redshift 
sources may be a viable way to gain further knowledge 
of cosmological parameters although there are still unre- 
solved issues concerning the modeling of the lenses. E.g. 



compact, the derived value of Hq is too low in comparison 
with the local measurements. 

In this note we investigate how to use multiply imaged 
high-redshift supernovae (SNe) to constrain the Hubble 
parameter, the mass and dark energy density of the uni- 
verse, f^M and fix, as well as the equation of state pa- 
rameter, Wq = which we assume is constant for z<,b. 
SNe are well suited for this technique because of the ex- 
pected high rate at high-redshifts and most importantly 
because of their well known lightcurves. In particular, the 
(rest-frame) optical lightcurves of core-collapse (CC) SNe 
show fast rise-times, typically about 1 week long. Thus, 
the time difference between images can be measured to 
better than one day's precision. 

Another possibility is the use of the UV shock- 
breakout. As has been seen in SN 1987A and modeled in 
Ensman fc Burrows (1992] ) the shock breakout can serve 



in a recent paper ( Kochaneck 2002 ) , five gravitational lens 



systems are analyzed and it is suggested that, unless the 
dark matter distributions in the lensing galaxies are rather 



as a time stamp with a precision of just minutes. The en- 
tire UV flash occurs over a period of minutes to several 
hours in the rest-frame of the SN depending on the nature 
of the progenitor. 

Wide field optical and NIR deep surveys, such as the 
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planned SNAP satellite (Pcrlmuttcr ct al. 2001), have the 
potential to discover ~ 10 B CC SNe. While these SNe are 
sometimes regarded as a "background" for the primary 
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cosmology program based on Type la SNe, we argue that 
lensed SNe of any kind may also provide useful informa- 
tion on cosmological parameters. 



In Holz (2001) the number of multiply imaged CC SNe 
up to z < 2 were estimated by simply scaling the Type la 
rate by a factor of 5. In this work we extend the considered 
redshift up to z — 5 using a SN rate calculation derived 
from the star formation history. Further, we take into ac- 
count the NIR wide field instrumentation in the current 
design of the SNAP satellite. 

Using simple toy-models, we investigate the accuracy 
of the strong lensing technique to improve our knowl- 
edge of cosmological parameters. Our observables are the 
source redshift, the redshift of the lensing galaxy, the 
image-separation, A8, the time-delay, At, and the flux- 
ratio, r. We derive the relation of our observables and 
cosmological parameters for different matter distributions 
in order to investigate the sensitivity to the choice of lens 
model. We speculate, that if the measurement is done us- 
ing a very large sample of lensed systems, useful bounds 
on the cosmological parameters can be found in spite of 
large uncertainties in the lens model. 

We have used the SNOC Monte-Carlo simulation pack- 
age ( poobar et al. 2002b ) to estimate the rate and measur- 
able quantities of multiple image SNe, e.g. the distribution 
of time-delays, image-separations and flux-ratios. We also 
simulate extinction by dust both in the host galaxy of the 
SN and in the foreground lensing galaxy. 

While about half of the known lensed systems exhibit 
more images (see e.g. Keeton et al. 199§| ), this work is lim- 
ited to spherically symmetric lensing systems producing 
only two or ring-like images. Systems with more lensed 
images are potentially very interesting as they provide 
more measurables which can be used to constrain the lens 
model. 

2. Modeling of galaxy halos 

In this work, we confine our study to a very simple class 
of spherically symmetric le ns models with project ed, two- 
dimensional lens potential ( Bchncider et al. 1992 ) 



ip(r) = kr a , 



(1) 



where a < 1 to assure that we get multiple images and 
that the surface density falls at large radii. The two limit- 
ing cases of this class of models are the Singular Isothermal 
Sphere (SIS) where a — 1 and the point-mass for which 
a = and ka = 1. A drawback with this choice of single- 
slope models is that we need steep halo models in order 
to produce multiple images. A slightly more general treat- 
ment would include double slope models like the isother- 
mal sphere with a core (ISC) or the NFW density profile 
(Navarro et al. 1997). Also, since many multiple image 
systems have more than two images, the ellipticity of the 
lensing galaxy should be included in the model. That is, 
since we assume a very limited class of lens models, we un- 
derestimate the systematic effects from the lens modeling. 
However, we also use a very limited set of observables. In 



four-image systems, we will have ten observables to use in 
the fit (four image positions, three flux-ratios and three 
time-delays) in comparison to the three observables used 
in this study. Still, there is no doubt that we underesti- 
mate the errors in the lensing potential while restricting 
our study to lens potentials of the form given by Eq. ([!]) . 
As a first step, we compare results for the two extreme 
cases of point-masses and SIS density profiles. 

2.1. Point- masses 

A light ray which passes by a point-mass M at a minimum 
distance £, is deflected by the "Einstein angle" 



AM 
T 



2R S 



(2) 



where Rs = 2M is the Schwarzschild radius. Using a char- 
acteristic length £o i n the lens plane given by 



2R s £^=p oDd , 



2R S 



Dds 
DdDs' 



(3) 



(4) 



where Dd s , T)d and D s are angular diameter distances be- 
tween deflector and source, deflector and observer, and 
source and observer respectively. The lens equation in di- 
mensionless form is 



y = x- 1/x, 

which has two solutions, 



Xi 



(5) 



(6) 



i.e., one image on each side of the lens. The image sepa- 
ration is given by 



A0 = /W?/ 2 + 4. 
The magnification is given by 
1 



which can be combined with Eq. ^| to give 
f 



Ml,2 = ±- 



VV +4 



±2 



y 



(7) 



(8) 



(9) 



The ratio r of the absolute values of the magnifications 
(flux-ratio) for the two images produced by a point-mass 
lens is given by, 







y/y 2 + 4 + y 






vV + 4 - y 



and 

y = r 1 ' i -r- 1 ' A . 



(10) 



(11) 
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Thus, we can write the image separation in terms of the 
flux-ratio as 



AO = /? (r 



i/4 , _-l/4 



)• 



(12) 



The time delay between images is in the general case 
given by 



DrlD. 



d^ds 



(l + z d ) 0(x^,y)-0(xW,y) , (13) 



(2) 



where <p(x,y) denotes the so called Fermat potential. In 
the case of a point-mass lens we have 



At = 2R s {l + z d )T(y), 
where 



ny) = ^vVv + 4 + m ^ — = — ■ 

2 y/y 2 + 4-y 



(14) 



(15) 



We can now express R$ and y in our observables to get 

(16) 



At = — (1 + z d ) 



2 D. 



ds 



i/V^ + 2 



2.2. Singular isothermal sphere (SIS) 

The density profile for the SIS halo, which is frequently 
used in gravitational lensing analysis, is given by 



v 2 1 
PSISO") = 2' 



(17) 



Here, v, the only free parameter of the model, is the line-of- 
sight velocity dispersion. The simplicity of the SIS density 
profile allows analytical solutions for the many quantities 



related to gravitational lensing (Turner et al. 1984). The 
deflection angle is computed to be 



2 lei ici 

— a - 



(18) 



where ao is defined as the magnitude of the deflection. 
That is, the magnitude of the deflection angle is indepen- 
dent of the impact parameter. Choosing 



£ =4tt 



v\ 2 D d D ds 



the lens equation can be written 

x 



y = x 



(19) 



(20) 



Considering y > 0, there are two images for any y < 1; at 
x = y + 1 and x = y — 1, i.e., on opposite sides of the lens 
center. The image separation is given by 



A6 = 2a, 



ds 



D 



and the magnification for an image at x is 



x-1 



(21) 



(22) 



From Eq. (|22j) we see that the ratio of the absolute values 
of magnifications of the images is 

1 



y 



1-1/ 

The time delay for the two images is 



(23) 



At 



4tt 



D„ 



0- + Zd)2y, 



(24) 



which can be expressed in terms of our observables as 



At 



A9 2 D d D s 
2 D ds 



(l + z d )- 



1 



1 



(25) 



We see that only the r-dependent part differs between the 
point and the SIS case. It is reassuring to see that they 
only differ by at most 6% up to a flux ratio of four. Thus, 
the derived value of the cosmological parameters are only 
weakly sensitive to the choice of point or SIS model. 

3. Monte- Carlo simulations 

Since we in principle only are able to follow infinitesi- 
mal light-beams in SNOC, we have to use some approx- 
imations when trying to get information on multiple im- 
age systems. The main approximation is that we assume 
that in cases of strong lensing, the effects from one close 
encounter is dominant, i.e., the one-lens approximation. 
With this simplification, we can use the information from 
the magnification to derive quantities for systems with fi- 
nite separations. In order to do this, we need to be able to 
derive analytical relations between the magnification and 
the image-separation and so forth. Here we show how this 
is done for the case of SIS lenses. 

First, we concentrate on primary images. Studying 
Eqs. ( pp| ) and j2^), we see that multiple imaging occurs 
whenever [i\ > 2 and that the magnification of the second 
image is given by 

M=Mi-2. (26) 

Using Eq. ( [HI) and the fact that y — l/(/xi — 1), we can 
write the time-delay for a SIS lens as 



At = 



47T 



(- 

V c 



DdDds 
D„ 



{l + z d )- 



1 



(27) 



Therefore, in order to compute the quantities of interest, 
we need to pick a lens redshift and velocity dispersion 
from some reasonable distributions for every case where 
/ii > 2. The distributions will in general be functions of 
the cosmology, the mass distribution of the lenses and the 
source redshift. The differential probability for multiple 
imaging is in the general case given by 

dP cx a{z d , z s )-^dM(l + z d ) 3 dV, (28) 

where o~(z d ,z s ) is the cross-section for multiple imaging 
and n is the comoving number density of lenses. For SIS 
lenses, the cross-section is given by 

*(z d ,z s )<x {-^ D 2 ds . (29) 
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bmce dV oc D 2 d -^dz d and n is independent of the lens 



redshift we can use 



dP(z d )^D 2 ds D 2 d (l + z d ) 3 ^-, 



(30) 



as our probability distribution for z d . The probability dis- 
tribution for v is given by 



. dn dM /v\ 4 , 
dPM oc — — — ; — - dv. 

w dM dv \c 



(31) 



dn = (f){y)dy 

<t>iv) = 



y a e~y, 



(32) 



with the mass-to-luminosity ratio (Peebles 1993, Eq. 3.39) 
M 



Using Eq. (||), we find that 



dn 



dM 

5 = a 



oc y s e~ v , 



ft 



1-/3 

Combining the Faber- Jackson relation 



(34) 
(35) 

(36) 



where y is defined in Eq. (|32|), with the mass-to- luminosity 
ratio, Eq. (|33|), we can relate the velocity dispersion to 
mass by 



v 



M 7(1-/3) 



(37) 



In this paper we use v* = 220 km s _1 , a — —1.07, ft = 
0.2 and 7 = 0.25. The mass normalization is calculated 
assuming that the entire mass of the resides in galaxy 
halos (Bergstrom et al. 2000). 

A higher value of u* would result in a larger number 
of wide separation lenses since the cross-section for mul- 
tiple imaging scales as a oc v 4 [see Eq. (29)] and the im- 
age separation as A9 oc v 2 [Eq. (pl|)]. Varying the galaxy 
mass distribution will also have an effect on the charac- 
teristics of the lensed events. However, since the mass and 
velocity dispersion are not variables in the fitting of cos- 
mological parameters (see Sect. FJ), any changes in these 
distributions will only have a marginal effect on the results 
obtained in this paper. 



700 



600 - 



500 



Z400 h 



Following Bergstrom et al. (2000| ) , we derive a galaxy mass J*300 - 
distribution, dn/dM , by co mbining the Schechter luminos- 
ity function ([Peebles 1993|, Eq. 5.129), 



200 - 



100 



• J— bond threshold 
O I— bond threshold 




25 25.5 26 26.5 27 27.5 28 28.5 29 29.5 30 
Magnitude threshold 



Fig. 1. Number of multiple image SNe in the simulated 
(33) sample vs magnitude threshold for the faintest image in I 
and J-band. 



4. Characteristics of lensed SN events 



We have simulated 1.1 • 10 6 CC SNe in the redshift range 
< z < 5 using a SIS model for the galaxy density profiles. 
The redshift dependence of the SN rates, the relative frac- 
tion among the various types of CC SNe (Ib/c, IIL, Iln, 
IIP, 87a-like), peak magnitudes and intrinsic dispersion 
were simulated following the prescriptions in Dahlen & 
|Fransson (1999TJ and |Dahlen fc Goobar (2002] ). Using SN 
rate predictions derived from the star formation history 
Sullivan et al. 200C), our simulated 



(Chugai et al. 1999 



sample corresponds approximately to the predicted num- 
ber of CC SN explosions (z < 5) in a period of 3-years in 
a 20 square degree field, i.e., ~ 5.1 SNe year -1 arcmin~ 2 . 
In the simulations we have assumed a flat cosmology with 
f^M = 0.3 and = 0.7, and a Hubble constant Ho = 65 
km s _1 Mpc _1 . We view these rates on the conservatively 
low side given the recent work by Lanzetta et al. (2002| ) 



who finds that the SFR plausibly increases monotonically 
with redshift through the highest redshifts observed. 

Out of the 1.1 • 10 6 simulated CC SNe, 2613 were mul- 
tiply lensed. Fig. [I] shows the number of detectable SNe 
as a function of the peak-brightness threshold in I and 
J-bands. The threshold refers to the faintest of the two 
images. 

In the following analysis, we consider two possible data 
samples (all magnitudes are in the Vega system): 

— case A: only SNe with an I-band peak magnitude of 
the faintest image of mi < 28.5 
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case B: SNe where either the I-band peak brightness 
or the J-band brightness falls below 28.5 mag (i.e. mi 
or m, 7 < 28.5) 



Thus, only 211 (case A) and 366 (case B) lensed CC 
SNe were considered for cosmology fits. Additional poten- 
tial SNe discovered in other bands are thus not considered 
in this study. For an instrument like SNAP (wavelength 
sensitivity 0.4-1.7 /im), revisiting each field in at least one 
optical or NIR band every 6 (or less) hours, we expect 
a very high detection efficiency. E.g., assuming that the 
SNe were discovered 0.5 magnitudes below peak bright- 
ness, SNAP would be able to make a 5a detection every 
four days in each of three (or more) filters: I, Z , J and 
possibly others. 

In the simulations we have considered the possibility 
that one or both images suffer extinction. We assumed dif- 
ferential extinction properties as in |Cardelli et al. (1989| ), 
with Ry = 3.1 and a mean- free path of 1 kpc for V-band 
photons in the host and lensing galaxy. Details of the sim- 
ulation procedure can be found in |Goobar et al. (2002a| .|b|) . 

In Figs. we show the distribution of observables 
for the multiple-image SNe; source redshifts, image flux- 
ratios, image separations and time-delay between the im- 
ages. The dashed lines show the subset of the data fulfilling 
case A SNe while the dotted curves correspond to case B. 

The salient features of the events are: images of the 
same source with average redshift z s ~ 2 and the lensing 
galaxy typically at z d ~ z s /2. The SN images are sepa- 
rated by ~ 0.5 arcseconds and a few weeks apart. With the 
imposed minimal brightness criteria, the images differ typ- 
ically by less than one astronomical magnitude. Clearly, 
a very clean signature to search for, especially with space 
observations free of atmospheric blurring. 



5. Fitting cosmological parameters 

Next, we estimate the target statistical uncertainty of a 
3 year mission studying a 20 square degree patch of the 
sky down to 28.5 peak I-band (and J-band) magnitudes. 
We perform a maximal likelihood analysis of the multi- 
ply lensed systems with h, f^M and Q\ as free parameters 
where h is the dimensionless Hubble parameter defined 
by H = h-100 km s -1 Mpc _1 . With = (h,n M ,^\) 
and f(r) being the r-dependent factors in Eqs. ( |l6| ) and 
( p5| ) for point mass lenses and SIS respectively we define 
the "experimental" (i? cxp ) and cosmology dependent (1Z) 
quantities: 



n{<j>) = 



and 



100 - 



80 - 



60 - 




40 - 



20 - 



2 2.5 3 
Redshift 



Fig. 2. Redshift distribution of lensed SNe up to z = 5, 
solid line. The dotted line shows SNe for which the SN 
system satisfies the case B criteria. The dashed lines shows 
case A. A SIS lens profile was used. The lensing galaxy is 
typically at a redshift ~ z s /2 

The maximum likelihood analysis performed over the 
simulated SNe with index i =1„N assumes Gaussian prob- 
ability density functions for the observables, 



i 



(38) 



where the standard deviation a contains the propa- 
gated uncertainties of the (independent) observables x = 
(At,A0,r, z d ): 



dR c 



4 

^ ^ <).<■ ; 



(39) 



D 



ds 



For simplicity, we have assumed that all SNe are measured 
with the same precision. 

While we will consider more pessimistic scenarios later, 
in this section we assume that the uncertainties o~ Xj of the 
four observables are: 



(40) 



The outstanding precision in the time delay between 
the SN images is consistent with the estimates of the ac- 
curacy expected for SNIa lightcurves studied with simu- 
lations by the SNAP collaboration. 

The image separation uncertainty corresponds to 0.1 
pixel of the proposed SNAP instrument. This estimate 



0~At 


= 0.05 days 




= 0.01" 


0~ r 


= 0.1 




= 0.001 
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Fig. 3. Flux-ratio of SN images. The dotted line shows 
SNe for which the SN system satisfies the case B criteria. 
The dashed line shows case A. A SIS lens profile was used. 



Fig. 4. Image separation in arcseconds. The dotted line 
shows SNe for which the SN system satisfies the case B 
criteria. The dashed line shows case A. A SIS lens profile 
was used. 



could even improve considering that all of the SN images 
along the lightcurve can be co-added t o get >40cr signals 
on the individually lensed SNe. E.g. in Anderson & King 
(2OOO0it is stated that with HST/WFPC one could reach 
0.02 pixel precision on reasonable bright stars. 

In estimating oy we have assumed that the differential 
extinction of the lensed SNe will not largely exceed what 
has been measured for a set of 23 gravitational lens galax- 
ies in Falco et al. (1999). The median extinction for those 
systems (z^l) was found to be E(B — V) ~ 0.05 mag. 

5.1. The f2fvi — plane 

While 1Z cx h , it varies only weakly with Qm and fix- 
This can be appreciated from the projected 68 % confi- 
dence level (CL) region from a 3-parameter fit (ilyi, ^A, h) 
to the simulated data-set B in Fig. ^. From now on we 
concentrate on the case B sample. The areas covered 
by the contour regions of case A scale approximately as 
^JNb/Na- It should be noted that the shape and inclina- 
tion of CL region is very different from what is found from 
the magnitude-redshift test ( Goobar fc Perlmutter 1995| ). 



We also show in dark the constraints in the il^ — 
plane if h is known (exactly). While little information can 
be gained from Qm or Q\ independently, the figures indi- 
cate that there is only a very limited range in 17m and fl\ 
for which the CL region is consistent with a flat universe, 
J7m + f^A = 1- Thus, from now on, we concentrate on the 
cases where the geometry is assumed to be known from, 
e.g., CMB anisotropy measurements. 



5.2. The — h plane 

With large statistics, such as the case in our simulations, 
it is feasible to constrain the possible values of J7m if a flat 
universe is assumed. Fig. [?| shows the 68 % CL-region of 
the ilu — h plane for case B. The dark shaded region shows 
the systematic effect introduced by fitting the SIS simu- 
lated data with the assumption of point mass lenses. The 
fit in the light shaded region (solid curve) of Fig. |t] yields 
h = 0.65±g go|; fl M = 0.30t^5. The uncertainties are 68 
% CL for a two parameter fit, i.e., 1.51 a. Introducing 
the bias due to the wrong lens model, we found the fitted 
central values to be h = 0.69 and £1m = 0.19. To further 
test potential bias in the cosmological parameters from 
systematic uncertainties in the lens model we have tested 
adding an ad-hoc offset with a different dependence on the 
image ratio: 



-^cxp — ^?cxp(l + £7" ), 



(41) 



where e = ±0.001 and ±0.005 and the considered un- 



certainties in the fit are as in Eq. 40 . The effect is shown as 
dotted and dashed curves in Fig. Note, however, that in 
these examples we are assuming that all lensing galaxies 
are wrongly modeled. Thus, in that sense, it should be a 
conservative estimate of the systematic uncertainties. 

It is interesting to note that if the Hubble parameter 
and the energy densities are known to good accuracy from 
other cosmological tests, it should be possible to put very 
useful constraints on the average lens properties. 



Goobar, Mortsell, Amamillah & Nugent: Cosmological parameters from lensed supernovae 



7 




60 80 100 120 140 160 180 200 
Time-delay (days) 



Fig. 5. Time delay between the SN images in days. The 
dotted line shows SNe for which the SN system satisfies 
the case B criteria. The dashed line shows case A. A SIS 
lens profile was used. 




Fig. 6. 68 % CL region in the — projection of a 
three parameter fit (f^M, ^A, h). The 366 two-image events 
fulfilling criteria B. The dark (green) region shows the 
68 % CL region that would result if h would be exactly 
known from independent measurements. The line shows 
J7m + ^A = 1, i.e. a flat universe. The diamond shows the 
value that was used in the simulation. 



5.3. The £1m — and h — wq planes 

Next, we investigate the potential of this method in setting 
limits in the Qm — wq parameter plane. Fig. ^| shows that 
upper limits on the equation of state parameter of dark 



0.75 



0.70 - 



A 0.65 - 



0.60 



0.55 



0.00 0.10 0.20 



0.30 



0.40 0.50 0.60 



Fig. 7. The light (yellow) shaded region bounded by the 
solid line shows 68 % CL region of J7m — h fit from 366 
two- image events where the peak brightness of the faintest 
image fulfills criteria B. A flat universe was assumed. The 
dark (green) shaded region shows the bias introduced if 
the fit is done (erroneously) assuming a point-mass lens 
model. The dotted and dashed curves show the effect of a 
systematic error in the lens model, according to Eq. [4l| 

energy^], u>o < —0.8 may be derived from the considered 
data sample. Meaningful limits on the possibility of wq < 
— 1 could be derived, especially if an independent estimate 
of Hq is used as a prior in the fit. 

Due to the different reshift distributions considered, 
the shape of the CL-regions differ from what is expected 
from the Hubble diagram of Type la SNe for the SNAP 
satellite ( Goliath et al. 2001 ). In particular, if the Hubble 
parameter is further constrained by some independent 
method, the estimates of wo from the strong lensing data 
becomes comparable in precision the limits that may be 
derived from Type la SNe, the dashed countour line in 
Fig. ||. The other 2D projection, the one onto the wq — h 
plane, is shown in Fig. ^. If the wrong halo model is used 
in the fitting procedure, i.e., a point-mass halo model in- 
stead of SIS, a 5 % bias is introduced in the estimate of 

Wq. 

6. Summary and conclusions 

Strongly gravitationally lensed SNe could be detected in 
large numbers in planned wide field, deep, SN search pro- 
grams, probably on the order of several hundred. In par- 
ticular, we argue in favor of large NIR imagers in space 
missions like SNAP. 

Lensed SNe are potentially interesting as they pro- 
vide independent measurements of cosmological parame- 
ters, mainly H , but also the energy density fractions and 
the equation of state of dark energy. The results are inde- 
pendent of, and would therefore complement the Type la 

1 Shown in the figure is the 2D projection of the 3-parameter 

fit (h, tt M , Wq). 
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fulfilling criteria B. A flat universe was assumed. The dark 
(green) region shows the smaller confidence region that 
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dent measurements. The dashed line shows the expected 
statistical uncertainty from a 3 year SNAP data sample 
of Type la SNe. 
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Fig. 9. 68 % CL region of h — wo fit from 366 two-image 
in case D. A flat universe was assumed. The dark 



events 

(green) region shows smaller confidence region that would 

result if £Im would be exactly known from independent 
measurements . 



program. At the faint limits considered in this note, sev- 
eral quasars per square arcminute are expected (Conti ct 
al. 199|). Thus, we expect that an instrument like SNAP 
would find several hundred multiply imaged QSOs, in ad- 
dition to the strongly lensed SNe. Thus, the statistical 
uncertainty could become smaller than what we have con- 
sidered here. 

While the systematic uncertainties remain a source of 
concern, we show that the simplest spherically symmetric 
models introduce moderate biases (ctq^^O.I ) , at least as 
long as multiple images with similar fluxes are considered. 



While projections of different SNe could in principle be 
interpreted as a lensed SN the two scenarios may be dis- 
tinguished. The signatures of CC SNe are unique. Crudely, 
the lightcurves are a product of the progenitor mass, the 
mass loss during its evolution off the main sequence, the 
amount of radioactive Ni synthesized during the explosion 
and the kinetic energy imparted to the ejecta. In addition, 
the environments the SNe explode in (the density and 
structure of their local interstellar medium) , often play a 
significant role in what we eventually see of the CC event. 
These differences have given rise to all the different classi- 
fications of these events we currently have; Type IIP, IIL, 
lb, Ic, Iln, etc. Given all this diversity it makes it quite 
easy to distinguish one event from another and to not con- 
fuse a lensing event with a coincident CC event along the 
same line of sight. 

In our simulations we find that extinction in the fore- 
ground galaxies does not severely affect the detectabil- 
ity of multiple lensed events nor the ability to derive the 
flux-ratio between the images. Clearly, multi-band obser- 
vations of the SNe will be important in order to correct 
for the different amounts of extinction of the images. At 
the same time, the data could provide important results on 
the dust properties of the foreground lensing galaxies, sim- 
ilar to the studies done with multiple imaged quasars at 
Zd^iA ( Falco et al. 1999 ). Unlike quasars, CC events have 
a very deliberate color evolution along their lightcurves. 
In general, from the moment of shock-breakout onwards, 
the atmospheres of CC SNe expand and become cooler 
and redder. This signature not only helps in the relative 
timing of these events, but also allows us to measure the 
differential extinction due to the varying amounts of dust 
along the lensed paths to the SN quite well. With 3 or 
more filters one could make a measurement of both the 
total extinction relative to the bluest event in addition to 
the ratio of the total to selective extinction due to differ- 
ences in the dust properties. Furthermore, one could en- 
hance this method by taking a spectrum of the SN (any of 
the lensed events would do) at a given epoch and through 
spectrum synthesis derive the true, unextinguished spec- 



tral energy distribution of the event [see e.g. (Mitchell ct 
al. 2002|; |Baron ct al. 2002|)]. 
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